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ABSTRACT 

A  simplified  analysis  is  performed  of  the  turbulent  mixing 
of  a  trace  chemical  species  with  the  bulk  chemical  species.  The 
general  relationship  between  the  mixing  on  a  molecular  scale  and 
the  concentration  fluctuation  intensity  of  the  trace  species  is  given. 
This  is  then  applied  to  a  turbulent  wake,  to  determine  the  axial  length 
that  it  takes  for  injested  gas  to  become  mixed.  The  results  are 
sensitive  to  the  assumed  turbulent  dissipation  rate;  giving  values 
of  the  mixing  length  from  the  order  of  to  two  orders  of  magnitude 
larger  than  the  wake  radius. 
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Introduction 

In  performing  calculations  of  chemical  reactions  in  turbulent  wakes, 

it  is  necessary  to  account  for  the  rate  at  which  the  gas  ingested  by  the 

turbulent  front  becomes  mixed  on  a  molecular  scale  with  the  gas  on  the 

interior  of  the  turbulent  wake.  One  proposed  method  uses  a  "mixing  lag 

length"  L,  which  is  deduced  on  the  basis  of  isotropic  turbulence  theory;^ 

another  method  uses  a  transport  equation  for  the  electron  density  fluctua- 

2 

tiens  to  determine  the  mixing  rate.  In  this  note,  it  is  shown  that  the  latter 
equation  caui  be  used  to  determine  the  mixing  length. 

* 

Analysis 

The  model  for  determining  this  length  is  a  one -dimensional  wake 
of  uniform  properties  of  species  A  which  ingests  a  small  amount  of  species 
B  at  a  distance  Xj  behind  the  body.  The  flow  field  exterior  of  the 


wake  is  also  assumed  to  consist  only  of  species  A. 

A  mixing  parameter  0  is  defined  as  follows; 


^  CA  Cb/CA  (1) 

t 

where  and  Cg  are  the  mass  concentrations  of  species  A  and  B 

respectively,  and  the  overbar  refers  to  a  time  average  at  a  particular  sta¬ 
tion  in  body-fixed  coordinates.  If  the  gas  is  ur.mixed,  then  either  CA  =  1, 

Cg  =  0,  or  CA  ~  0,  Cg  =1,  so  that  Cg  =  0;  that  is  0  =  0.  If  species  B 
is  intimately  mixed  with  A,  then  Cg  (t)/CA  (t)  =  constant  =  ^g/C^»  but  s*nce 
then  (t)  =  CA,  0=1.  Thus  0  is  a  measurement  of  mixing  on  a 
molecular  scale.  Also,  since  there  are  only  two  species  in  this  model 


cA(t)  +  Cg  (t)  =rr,  +£b=  i 


(2) 


If  we  decompose  each  concentration  into  a  mean  value  and  a  fluctuating 

t 

component  (  )  ,  from  eq.  (2) 


i  i 

C.  +  C_  =  0 
A  B 


(3) 


Combination  of  eqs.  (1,  2,  3)  yields; 


0=1—  Cg  ^ 


A  ~B 


(4) 


If  B  is  a  trace  species,  I>  and  eq.  (4)  can  be  rewritten  as; 


0(x)  =  1  - 


C„  (x)  TJg  (xj) 


B 


V  <*i> 


Cb2(x,|  CB  w  rB  (Xj) 


(5) 


where  x  is  the  distance  downstream  and  is  the  point  of  ingestion  of 
species  B  which  is  initially  unmixed  at  the  point  of  ingestion;  hence 
<fi  (Xj)  =  0  so  that  the  last  factoi  in  eq.  (5)  is  equal  to  unity,  e.g. 


cj2  (Xj)  ="CB  (x^  (6) 

The  second  factor  in  eq.  (5)  is  determined  only  by  dilution.  Thus  the 

problem  is  reduced  to  that  of  determining  ihe  decay  of  mean  square  fluctua- 

2 

tions  of  species  B,  for  which  Lin's  equation  (14)  can  be  used,  e.g.: 


(d/d  i)  ( r ™  u 


ufd  (rfm)/d£  =  -  Rr 


.1/3  k  2/3 

e 


m 


'B 


(7) 


where  m  =  1  for  plane  wakes  and  2  for  axially  symmetric  wakes;  r£  (x) 
is  the  mean  distance  from  the  axis  to  the  wake  turbulent  front;  £  =  x/^n’ 
Rn  is  the  wake  drag  or  momentum  defect  radius,  u  is  the  mean  wake 
speed,  e  is  the  rate  of  dissipation  of  turbulent  kinetic  energy  and  kg  is 
the  turbulent  energy- containing  wave  number.  The  dissipation  rate  given 
in  eq.  (7)  can  be  deduced  from  eqs.  (7,  11,  20,  and  23)  of  ref.  3,  for  a 
Schmidt  Sc  number  of  unity.  This  form  is  also  applicable  for  Sc  <  1,  if 
multiplied  by  ^  (3  -  S  ^). 


— HT  .  -*• 

From  eq.  (6),  Cg  is  of  the  order  of  CR  .  Since  CR  «  1, 


B 


2  *2 

«  Cn  and  the  second  term  in  eq.  (7)  can  be  neglected.  In  addition, 
±3  D 

one  may  approximate  U£  as  u.  Then  eq.  (7)  can  t>e  integrated  to: 

>: 

L  [r£m  CB  ?7'rfm  (x£)  (x£)J=  -  Rr  J  u  fc1'  '  k^3  d£  (8) 


-3- 


Now,  from  conservation  of  specie^  B,  r^‘  (x^)  (\j),  and 

use  of  eq.  (6)  m  eq.  (8)  yields: 


X  (x)  =  1  -  exp 


; -  Rr  /  u '*  ei/3  ke2/3  de ) 


(9) 


We  will  next  apply  eq.  (9)  to  throe  cases:  (1)  self-preserving  wake;  (2)  the 

2 

turbulence  parameters  *  alc  ilated  by  Lin  ;  and  (3)  the  turbulence  parameters 

4 

deduced  by  Schapker  from  turbulent  front  measurements. 

^elf  Preserving  Wake 

The  rate  of  destruction  of  turbulent  kinetic  energy  for  a  vake  can  be 
estimated  from  the  generalized  form  of  eq.  (8.4)  of  io£.  5  for  an  incompressible 
wake :  *, 

OO  M 

r  f  r  2  —*.*1 

m 


/ 


€  dr 


m 


1  d 

7  ^ 


r  r  _  ^  -r{\ 

I  u  (  Au  )  +  q  J  dr 


(10) 


o 


where  A  u  is  the  mean  velocity  defect  and  ^  q  2  is  the  turbulent  kinetic 
energy.  For  a  self  preserving  wake,  q'^  -  (Au)  .  Also,  from  conserva¬ 
tion  of  mass, 


m 


(10a) 


where  U  fs  the  free  stream  velocity.  Thus, 
00 


m  - 1 


€  (x)  ~  U  (R  /rjm  x 
'  oo  '  n  f 


(H) 


Since  the  self  preserving  wake  is  also  geometrically  similai,  the  energy- 
cont  ining  wave  number  is  xversely  proportional  to  the  wake  radius, 


-4- 


given  by 


1 

m  +  1 


(rf/Rn)  ~  %  (12) 

f 

and  the  mean  wake  velocity  is  taken  as  approximately  U^.  Use  of 
eqs .  (11,  12)  in  (9)  yields: 

<t>  =  i  -  (IjA)0  (13) 

where  c  is  the  constant  representing  the  numerical  constants  of  propor¬ 
tionality  in  eqs,  (11,  12)  and  between  kg  and  l/r^.  The  important  conclusion 
from  (13)  is  that  th  :  distance  required  to  achieve  a  given  amount  of  mixing 
depends  linearly  oa  the  distance  at  which  the  species  was  initially 
ingested:  that  is,  the  further  downstream  the  ingestion,  the  lc.»  ger  the 
length  to  become  well -mixed.  From  eq.  (13)  the  length  A£  =  £  -  Sj»  to 
achieve  a  certain  degree  of  mixing  is 

* 

-l/c 

'  AS/S,  =  (!-*>  -1  <14) 

1  I 

From  eq.  (14),  A|  —  °°  as  $  —  1  but  for  finite  values  of  (j>  and  c  ^/n(l-  <f>) «  1 

AS  *  -  c_1  Sj  U  (1-0)  (15) 

In  general,  the  integrand  of  eq.  (9)  can  be  represented  as: 

a  t;’1  £‘/3k  z^-  ci'n  (i6) 

a  n  e 


-5- 


and  eq.  (9)  can  be  integrated  as  follows: 


0=1-  exp 


f  -1  l~n  Un  1 

I-  =  1>  -  e  )J 


(17) 


For  either  c  large  or  0  small,  eq.  (17)  becomes; 


r;  -  c 


^l)  ^ 


(1-0) 


(18) 


which  is  the  general;.,  ation  of  eq.  (15).  In  order  for  the  mixing  length  A  £ 

to  be  a  constant  multiple  of  the  distance  from  the  axis  to  the  turbulent  front 

from  eq.  (12)  n  =  (m  +  1)  S  e.g.  1/2  for  plane  wakes  and  l/3  for  axially 

4 

symmetric  wakes.  This  condition  is  almost  fulfilled  for  Schapker's  data 
3 

for  £  <  6  x  10  ,  where  n  =  0.43,  as  shown  below. 

Enuation  (18)  can  also  be  written  as 


AC  =  -in  (1  -  0)[^Rn  u  e1/3  ke2//3j 


-  1 


e, 


(19) 


which  could  have  been  obtained  directly  from  eq.  (9)  by  taking  the  integrand 
constant  at  its  value  at  If  one  defines  a  characteristic  decay  time  T  as 


-i/3  -2/3 

7  «  €  k 


(20) 


and  a  non-dimensional  decay  length  A  =  U^t/R^,  then  eq.  (19)  becomes 


A£  =  -  A  /n  (1  -  0) 


(21) 


-  o  - 


which  yields  the  expected  exponential  decay  of  <f>.  The  use  of  a  definite 
mixing  lag  is  obviously  somewhat  arbitrary,  since  A  £  depends  on  the 
criteria  for  <j>.  To  compare  various  theories,  an  arbitrary  value  of  <j> 
equal  to  l/Z  has  been  chosen. 

Hypersonic  Turbulence  Parameters 

2  _  1 

Lin  assumed  k  ss  r  and  also  calculated  e.  A  good  match  to 
e  n  ° 

-  1 

those  results  for  £  >  10^  is  c  -  0.465,  n  =  2/3,  so  that  for  <f>  - 


A  4 


1.49 


(22) 


which  is  shown  in  Fig.  1  in  comparison  to  the  wake  radius.  It  is  seen  that 
the  mixing  is  between  one  and  two  orders  of  magnitude  larger  than  the 
mean  radius  of  the  turbulent  front,  which  is  primarily  a  consequence  of 

4 

the  low  values  of  f.  Schapker's  results  can  be  represented  as  follows 


--1  1/3  2/3  5/6  ,  0.4  -0.25  -0.43  . 

Rn  u  €  7  ke  /‘5=  86.5(3  X  (Tf/T„)  Redoo  £  ;  £  <  2000 


o  0  4s  /a  —0.25  _  ] 

=  6.75x  103  (Tf/Too)  (3  Re,  £  ;£> 


(23) 


2000 


where  (3  =  Uf/Uoo  ~  1  and  T  is  temperature.  Neglecting  the  factor  Tf/Too  , 
the  mixing  length  for  <j>  =  is  then  given  by: 


.2  0.25  0.43 

A£  =  1.15x10  (|n  2)  Redoo  ^ 


-4  »  0.25 

=  1.48  x  10  (in  2)  Red  £j 


e2<  2  x  io3 

>  2  X  103 


(24) 
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ING  LAG  ,A£ 


DIMENSIONLESS  DISTANCE  BEHIND  BODY, 

Fig.  1  Present  mixing  length  calculations  based  on  Lin' s^  and 

Schapker' turbulence  parameters  for  Re  ,  =  10^,  compared 

U  oo 

with  the  wake  radius  and  previous  estimates.  1 


A949? 


RADIUS 


and  is  also  shown  in  Fig.  1.  These  values  result  in  a  short  mixing  length, 

approximately  equal  to  the  wake  radius  within  about  a  factor  of  2.  These 

short  mixing  lengths  are  caused  primarily  by  the  large  values  of  both  e  and 

4 

k^  deduced  by  Schapker. 

Proudian  and  Feldman's  analysis^  results  in: 

A  =  §1  ^  -  exp  -  (p/p^)^3  n_1J  (25) 

1 

7 

where  p  is  mass  density  and  n  is  the  exponent  oi  the  defect  velocity  decay. 
For  a  self  similar  wake  n  ~  4/3,  and  for  p  =  pB,  <j>  = 

Ae  =  0.  37  il  (26) 

also  shown  in  Fig.  1.  Here  it  is  seen  that  this  mixing  length  is  much 
larger  than  any  of  the  others,  with  a  considerably  larger  slope.  The  values 
recommended  in  ref.  1  are  also  shown  in  Fig.  1,  adjusted  for  <f>  ~  ^ , 
which  has  an  even  larger  average  slope.  (For  £  >  600,  a  constant  mixing 
time  of  2  msec  .s  recommended.  ) 

Discussion 

The  mixing  length,  is  deduced  from  various  theoretical  calculations 
and  measurements  are  obviously  in  wide  disagreement  by  as  much  as  two 
orders  of  magnitude.  Unfortunately,  direct  measurements  are  not  avail¬ 
able,  but  calculations  have  been  performed  cf  the  effect  of  various  mixing 
lengths  on  electron  densities  in  wakes.  ^  The  calculations  shown  in.  Fig.  10 
of  ref.  6  indicate  that  the  results  are  insensitive  to  the  mixing  length  for 
£  <  2500,  but  becomes  very  sensitive  beyond  that. 


-9- 


Furthermore,  Figs.  3  and  10  of  ref.  6  which  compares  the  calcu- 
7 

lations  with  experiments  indicate  that  the  mixing  length  is  quite  short, 
e.  g.  A  i/i'f  S  10  for  x/Rj^  <  4000.  These  valves  are  less  than  those 
predicted  by  Lin  by  only  a  factor  of  about  3  at  x/R^  =  4000.  The  small 
values  predicted  by  Schapker  may  be  implausible.  Additional  experimental 
evidence  would  be  desirable  to  resolve  these  differences. 
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